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We present the different consistency relations that can be seen as variations of the well known 
Suyama-Yamaguchi (SY) consistency relation tnl (|/nl)^, the latter involving the levels of non- 
gaussianity /nl and tnl in the primordial curvature perturbation It has been (implicitly) claimed 
that the following variation: rNL(ki,k3) ^ (|) /NL(ki)/NL(k3), which we call "the fourth and the 
best-defined variety", in the collapsed (for tnl) and squeezed (for /nl) limits is always satisfied 
independently of any physics; however, the proof depends sensitively on the assumption of scale- 
invariance (expressing this way the fourth variety of the SY consistency relation as tnl ^ (|/nl)^) 
which only applies for cosmological models involving Lorentz-invariant scalar fields (at least at 
tree level), leaving room for a strong violation of this variety of the consistency relation when 
non-trivial degrees of freedom, for instance vector fields, are in charge of the generation of the 
primordial curvature perturbation. With this in mind as a motivation, we explicitly state, in the 
first part of this work, under which conditions the SY consistency relation has been claimed to 
hold in its different varieties (implicitly) presented in the literature since its inception back in 2008; 
as a result, we show for the first time that the variety TNL(ki,ki ) ^ (|/NL(ki))', which we call 
"the fifth variety" , is always satisfied even when there is strong scale-dependence and high levels of 
statistical anisotropy as long as statistical homogeneity holds: thus, an observed violation of this 
specific variety would prevent the comparison between theory and observation, shaking this way 
the foundations of cosmology as a science. In the second part, we concern about the existence of 
non-trivial degrees of freedom, concretely vector fields for which the levels of non-gaussianity have 
been calculated for very few models; among them, and by making use of the 5N formalism at tree 
level, we study a class of models that includes the vector curvaton scenario, vector infiation, and the 
hybrid inflation with coupled vector and scalar "waterfall field" where ( is generated at the end of 
inflation, finding that the fourth variety of the SY consistency relation is indeed hardly violated for 
some specific wavevector configurations while the fifth variety continues to be well satisfied. Finally, 
as a byproduct of our investigation, we draw attention to a quite recently demonstrated variety of 
the SY consistency relation: t^l (|/nl)^i in scenarios where scalar and vector fields contribute to 
the generation of the primordial curvature perturbation; this variety of the SY consistency relation is 
satisfied although the isotropic pieces of the non-gaussianity parameters receive contributions from 
the vector fields. We discuss further implications for observational cosmology. 
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I. INTRODUCTION: THE FIRST VARIETY OF THE SY CONSISTENCY RELATION 

Modern cosmology not only studies the background dynamics of any inflationary model of the Universe, but also 
the perturbation dynamics via the connected n-point correlators (C(ki)^(k2)...C(k„))c of the primordial curvature 
perturbation C [T|. During some time, the two-point correlator of ^ was enough for the purpose of comparing theoretical 
predictions with observation, but, after satellite missions reached amazing accuracy levels , it has been necessary to 
work out the connected three- and four-point correlators of ^. Non- vanishing connected three- or four-point correlators 
of C imply that the probability distribution function of the primordial curvature perturbation is non-gaussian [4j , that 
being the reason why the functions /NL(ki, k2, ka), and TNL(ki, k2, ks, k4) and gNL (ki , k2 , ka , k4), that parameterize 
the connected three- and four-point correlators of ^ respectively are called the levels of non-gaussianity. 

A few years ago, Suyama and Yamaguchi [5^ showed that /nl and tnl satisfy the following consistency relation: 



as long as the following conditions are satisfied: 

• Condition 1: The calculation of /nl £^nd tnl is performed at tree level in the diagrammatic approach of the SN 
formalism [H [B] . 

• Condition 2: The inflationary dynamics is driven by any number of slowly-rolling scalar fields. 

• Condition 3: The fields involved are gaussian. 

• Condition 4- The field perturbations are scale-invariant. 

Although /nl and tnl are not directly comparable in the general case because they are functions of the wavevectors, 
the latter expression is valid since, under the conditions previously stated, /nl and tnl are scale-invariant [T. We 
will call Eq. ([T]) the first variety of the Suyama- Yamaguchi (SY) consistency relation. 

Because of the role of the SY consistency relation of Eq. ([T]) at ruling out classes of cosmological inflationary models 
in the event that the inequality is observationally violated, it is reasonable to wonder if it is possible to relax some 
or all the conditions previously stated while preserving the SY consistency relation. Suyama, Takahashi, Yamaguchi, 
and Yokoyama [8 , and Sugiyama, Komatsu, and Futamase [3, made efforts in this direction, trying to take into 
account the one- loop corrections to the three- and four-point correlation functions (relaxing condition 1), but adding 
a new condition: 

• Condition 5: /NL(ki, k2, ka) is evaluated in the squeezed limit (ki — )■ 0) while TNL(ki, k2, ka, k4) is evaluated 
in the collapsed limit (ki -I- k2 — 0). 

This effort was consolidated by Sugiyama [10 who demonstrated that the SY consistency relation is valid if conditions 
2-5 are met and condition 1 is replaced by the following one: 

• Condition 1: The calculation of /nl and tnl is performed at all loop corrections in the diagrammatic approach 
of the 5N formalism [H |B] . 

However, in an elegant work by Smith, LoVerde, and Zaldarriaga [TT, it was shown that the first variety of the SY 
consistency relation is valid if condition 5 is met while conditions 1-4 are replaced by the following ones: 

• Condition 1: The calculation of /nl and tnl is performed non-perturbatively. 

• Condition 2: The inflationary dynamics is arbitrary. 

• Condition 2a (actually, a consequence of condition 2): The fields involved (if any) can be non-gaussian. 

• Condition 3: Statistical homogeneity of C, is preserved. 

• Condition 4- Statistical isotropy of is preserved. 

The statistical homogeneity of states that the connected n-point correlators in the configuration space are in- 
variant under space translations [1]; this implies that, in the momentum space, the connected n-point correlators are 
proportional to a Dirac delta function: 




(1) 



(C(ki)C(k2)...C(k„)), (X (53(ki + k2 + ... + k„) 



(2) 
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The above implies that the wavevector configuration is such that ki, k2, and k„ form a warped polygon. The 
statistical homogeneity is fundamental in cosmology since it allows us to compare theoretical predictions with obser- 
vations (via the ergodic theorem, see Refs. |31[T2]). In contrast, the statistical isotropy of ( states that the connected 
n-point correlators in the configuration space are invariant under space rotations [T; this implies that, in the momen- 
tum space, and once statistical homogeneity has been imposed, each connected n-point correlator is proportional to 
a function Mi^ called the n — 1-spectrum which is invariant under rotations in the momentum space: 

(C(ki)C(k2)...C(k„))e = (27r)3<53(ki + k2 + ... + k„)M<;(ki, k2, k„) , (3) 

where 

Mc(ki,k2,...,k„) =Mc(ki,k2,...,k„), (4) 

being hi = TZ k-i with TZ being the rotation operator such that — TZ in the configuration space, i.e. the n — 1- 
spectrum is the same no matter the orientation of the warped polygon. The statistical isotropy is the usual assumption 
in cosmology since the inflationary dynamics is usually driven by scalar fields only, which, due to its nature, do not 
exhibit preferred directions. 

Some time after the work by Smith, LoVerde, and Zaldarriaga, and by employing an alternative technique, Assassi, 
Baumann, and Green obtained the same result as the former ones under the same assumptions. However, both 
proofs rely on this additional assumption: 

• Condition 6: /nl and tnl are scale-invariant, 

as clearly stated in a later work by Kehagias and Riotto [1^ . The latter authors come back to the original conditions 
1 and 5, except for condition 2 which is reinforced and conditions 3 and 4 which are relaxed: 

• Condition 2: The infiationary dynamics is driven by any number of slowly-rolling Lorentz-invariant scalar fields. 

• Condition 3: The fields involved can be non-gaussian. 

• Condition 4- The field perturbations, in principle, might not be scale-invariant. 

By paying attention to the symmetries present during the de Sitter epoch, and employing the operator product 
expansion technique, Kehagias and Riotto show that the connected two-, three-, and four-point correlators of the field 
perturbations are scale-invariant in the squeezed (for the three-point correlators) and collapsed (for the four-point 
correlators) limits. Thus, the first variety of the SY consistency relation, under the conditions just stated above, is 
shown to be valid but just at tree level in the diagrammatic approach to the SN formalism. Kehagias and Riotto 
assume that their proof is valid, even including loop corrections, for a wide range of models since loop corrections seem 
to become important just for rather marginal cases [TS] (see however [IHIIIZ]); nevertheless, their proof at all loop 
corrections, employing the results by Sugiyama [lU], is still incomplete as long as the connected n-point correlators 
in the field perturbations, with n > 4, are not shown to be scale-invariant in the required limits. 

The proofs by Smith, LoVerde, and Zaldarriaga [TP, and by Assassi, Baumann, and Green [13], are by far the 
most general and lead us to wonder whether one or some of the conditions these proofs are based on are unnecessary; 
they also lead us to state, assuming that the proof in 14J is valid at all orders, that a violation of the first variety of 
the SY consistency relation would imply that new non-trivial degrees of freedom play a role during inflation. Among 
the conditions that could be unnecessary, we think that the statistical isotropy condition, and the scale-invariance 
condition for /nl and tnl; are the ones to investigate since the statistical homogeneity condition is essential to compare 
theory with observation [U [T^] and since the other three conditions are absolutely general. Among the possible new 
non-trivial degrees of freedom, the vector flelds are particularly interesting since they are suitable candidates to 
explain the apparent violation of statistical isotropy observed in the 5-year data from the WMAP satellite [18l [19] 
(see however |20|). given that they define inherently a preferred direction for the expansion, for the distribution of 
primordial fluctuations, or for both [1] HJ-dl]. The vector fields, the statistical anisotropy, and the scale-dependence 
of /nl and tnl are indeed related to each other since, even in the simplest case when the expansion is isotropic, vector 
flelds break the statistical isotropy in C (33j, and can make /nl and tnl strongly scale-dependent at tree level [M] 
and even more at higher levels in the perturbative expansion [3] [51]-|39|. 

II. RELAXING THE STATISTICAL ISOTROPY CONDITION: THE SECOND AND THIRD 
VARIETIES OF THE SY CONSISTENCY RELATION 

In this section, we are going to reproduce the proof by Smith, LoVerde, and Zaldarriaga [TT], step by step, but this 
time relaxing the statistical isotropy condition. As discussed around Eq. Q, after imposing statistical homogeneity. 
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the connected n-point correlator of ( is written in terms of a function M^(ki,k2, ...,k„), called the n — l-spectrum, 
which, if statistical isotropy is imposed, becomes invariant under spatial rotations of its arguments; this in turn 
implies that the spectrum P^(ki,k2) and bispectrum i?j(ki, k2, ka) depend in this case only on the magnitudes 
of the wavevectors involved (i.e., P|^(ki,k2) — P(;{k) and _B^(ki, k2, ka) = _B^(fci, fc2, ^3)) while the trispectrum 
T^(ki, k2, ks, k4) and so on still depend on the whole wavevectors. When relaxing the statistical isotropy condition, 
we must take into account that P(;(ki,k2) and P^(ki,k2,k3) do depend on the whole wavevectors and so does 
/NL(ki, k2, ks) (as well as rNL(ki, k2, ks, k4) and (7NL(ki, k2, k3, k4)); however, by virtue of the reality condition on 

the relation P|j(ki,k2) = P^(— ki,— k2) is met jH [31j. From now on, by virtue of the statistical homogeneity 
condition, we will call P|^(ki,k2) just as P|^(ki). 

The trispectrum of ri;;(ki, k2, ks, k4), is parameterized, by definition, in terms of the spectrum of i.e. Pf(k), 
and the levels of non-gaussianity TNL(ki, k2, ks, k4) and 5NL(ki, k2, k3, k4) in the following way [ICTl HT] : 

T<;(ki,k2,k3,k4) = rNL(ki,k2,k3,k4) [Pc(ki)P(;(k3)P(;(ki +k2) + 11 pcrmutations] 
54 

+ ^ffAfL(ki,k2,k3,k4) [P(;(ki)P(;(k2)Pc(k3) +3 permutations] . (5) 

This relation is greatly simplified in the collapsed limit ki +k2 — t- (which indeed implies ks +k4 — > and, therefore, 
ki — >• — k2 and k3 —^i) and by making use of the reality condition of ( and the fact that P^(k) = ^P^(k) where 
P^(k) is almost scale-invariant according to observations: 

hm rc(ki,k2,k3,k4)=4rNL(ki,k3)Pc(ki)Pc(k3)Pc(ki+k2). (6) 
ki+k2-i-0 

In turn, the bispectrum of ^, i?^(ki, k2, k3), is parameterized, by definition, in terms of the spectrum of C, i.e. P^(k), 
and the level of non-gaussianity /NL(ki, k2, k3) in the following way [42 l I43 j : 
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P(;(ki,k2,k3) = -/NL(ki,k2,k3) [P^{ki)Pt^{k2) + cyclic permutations] . (7) 



The above relation is simplified in the squeezed limit ki — >■ (which indeed implies k2 — > — k3) and by making use 

271 

fc3 



again of the reality condition of C and the fact that Pi;(k) = ^Pi^(k): 



lim Bc(ki,k2,k3) - ^/NL(k2)Pc(ki)Pc(k2) . (8) 
We will now define an auxiliary field called -P(k) as 

p(k)-^/ ^^k- C(kOC(k - k') 

where 65 is a narrow band of wavevectors which are very near some kg and is its volume: 



The next step is to calculate the correlator (P*(ki)P(ki)): 

1 r d^k.d^k^ (C(-ki)c(-kL + ki)c(k3)C(kL - k3)) 

L^ebs (27r)6 Pc(ki)Pc(k3) 



(P*(ki)P(kO)^i^ / 



V2 / ^o\s\^^^p(v^ (C(-ki)C(-kL + ki)C(k3)C(kL - k3)), 

JkiMsebs (27r)6 P^{ki)P^{k3) I 



s 

+ (C(-ki)C(-ki + ki))(C(k3)C(kL - k3)) + (C(-ki)C(k3))(C(-kL + ki)C(ki - k3)) 

+ (C(-ki)C(kz. - k3))(C(-ki + ki)C(k3))] , (11) 

where we have used the Wick's theorem [T21I31]. Thus, by employing the statistical homogeneity condition and the 
reality condition of (, (P* (ki)P(kL)) is finally written as 

(P*(k,)P(k,)) = i, / -^,3(o) ^c(-k.-k. +k,k^^ 

JkiMebs [^^) ^C(kijn(k3j 
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The collapsed limit for the trispectrum of C in the above expression is obtained when — ki — k/, + ki — >■ 0, i.e. when 
ki — > 0. Thus, by employing the collapsed limit for the trispectrum of ^ in Eq. ([6]), and taking the limit k/, — >■ in 
the above expression, we find 



We will now calculate the correlator (C*(kL)P(ki)): 

1 



^^rNL(ki,k3) + - 



(C(k,.)P(ki)) = 



Vs 



k2e6s 



d^k2 (C(-kL)C(k2)C(kL-k2)), 
(27r)3 P^(k2) 



(13) 



(14) 



The squeezed limit for the bispectrum of C in the above expression is obtained when k^ — >■ 0. Thus, by employing 
the squeezed limit for the bispectrum of ( in Eq. (Isl), and taking the limit k^ ^ in the above expression, we find 



lim 



(C(kz.)P(k^)) 
krA'o (27r)353(o) 

Finally, let's construct the covariance matrix 



/NL(k2) . 



(r(ki)c(ki)) (c(ki)p(ki)) 

(27r)353(o) V(C(ki)P*(ki)) (P*(ki)P(ki)) 



(15) 



(16) 



whose determinant is positive; indeed, such a condition, in the k^ — )■ limit, and by introducing Eqs. (13) and (15 1, 
reduces to 



1 

1/2 , 



(27r)6 



rNL(ki,k3) > 



6 1 
5T^ 



r/i 



nH 



1 



2P^{kL)Vs 



(17) 



where the last term in the right hand side of the equation goes to zero in the k^ — > limit. Thus, we arrive to the 
actual expression obtained by Smith, LoVerde, and Zaldarriaga |llj . before assuming scale invariance, but this time 
relaxing the statistical isotropy condition: 



(2^)6 ^NL(ki,k3) ^ 
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d^k. 



5 ./k.efc. (2^)3 



/NL(k2 



(18) 



When statistical isotropy is required, the argument of /nl in the above inequality must be replaced so that it becomes 
the wavenumber ^2 instead of the wavevector k2. We will call Eq. (18) the second variety of the SY consistency 
relation. 

It is not our intention to reproduce in this paper the proof by Assassi, Baumann, and Green '13], but we can say 
that, relaxing the statistical isotropy condition, and before assuming scale invariance, their proof leads actually to the 
consistency relation 



^^^^rNL(ki,k3)Pc(ki)Pc(k3) ^ 



fk2 

{2ny 



/NL(k2)Pc(k2 



(19) 



where the integrals are performed over the entire momentum space. When statistical isotropy is required, the argu- 
ments of Pq and /nl in the above inequality must be replaced so that they become the respective wavenumbers fc„ 
instead of the respective wavevectors k„, fixing in this way a small mistake in Eq. (2.13) of [13]. We will call Eq. 
( 19 ) the third variety of the SY consistency relation. 



III. THE SCALE INVARIANCE ARGUMENT: THE FOURTH AND FIFTH VARIETIES OF THE SY 

CONSISTENCY RELATION 

The previous section showed what the actual consistency relations are, between tnl and /nl, before assuming scale 
invariance. If scale invariance is guaranteed, as argued by Kehagias and Riotto |T3] in the case that the inflationary 
dynamics is driven by Lorentz-invariant scalar fields only, the second and the third varieties of the SY consistency 
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relation, Eqs. (181 and ([19]), reduce to the first variety in Eq. ([T]). However, why should it be like that in the general 
case? Could it be that, although the second and the third varieties are satisfied, the first variety is violated for those 
cases where there is no scale invariance?, i.e. is it possible to have the following relation?: 



rNL(ki,k3) < 



./NL(ki)/NL(k3) . 



(20) 



It is worth noticing that the previous expression corresponds to a violation of what we will call the fourth and 
best-defined variety of the SY consistency relation: 




.fNL(ki)/NL(k3) 



(21) 



which is just a direct generalization of the first variety, Eq. ([ij, 
easily inspired from the second and third varieties, Eqs. ([Tsl and ([191 



when there is no scale-invariance and whose form is 
Well, it is clear that, if the proof by Kehagias 



and Riotto [13] is valid to all orders, the violation of the fourth variety of the SY consistency relation would imply the 
presence of new non-trivial degrees of freedom (sec for example Refs. [T31I1S]); however, the opposite is not necessarily 
true, and that is the motivation of our mission in the next section: to show that, when involving vector fields, there 
exist models where the inequality in Eq. (21 ) is indeed hardly violated for some wavevector configurations. Meanwhile, 



what we can say is that the best-defined variety of the SY consistency relation is satisfied when the arguments of tnl 
and /nl are all the same even if the scale invariance is not guaranteed: 



(22) 




this can be proven by observing that, in Eq 
but making the arguments of tnl and /nl all the same 
relation. 



(18), bs can be made as small as we want, getting rid of the integrals. 
We will call Eq. 



( 22 ) the fifth variety of the SY consistency 
The discussions in Sections [T] and [Tl] of this paper show that statistical homogeneity is actually the only 
non-general condition for the validity of the fifth variety of the SY consistency relation when rNL(ki, k2, ka, k4) and 
/NL(ki, k2, ka) are evaluated in the collapsed and squeezed limits respectively; an observed violation of such a variety 
would lead to the conclusion that statistical homogeneity does not hold and that, therefore, theory and observation 
cannot be compared [4] leading to a strong shaking of the foundations of cosmology as a science [46] . 



IV. A CLASS OF INFLATIONARY MODELS INVOLVING VECTOR FIELDS: VIOLATING THE 
FOURTH VARIETY OF THE SY CONSISTENCY RELATION 



In this section we show a worked example in which we get a violation of the fourth variety of the SY consistency 
relation in Eq. (21 ) at least at tree level. We begin by describing a class of infiationary models which includes several 



popular models of inflation in the presence of vector fields; this will be done using a generic parametrization. 



A. A class of inflationary models involving vector fields in view of the SN formalism 



The class of models that we shall consider in this section can be parametrized using its first three correlation 
functions. We use the SN formalism [31 [Z] to do our calculations at tree level for all the correlators and we will assume 
isotropic expansion. In the following, we shall use the notation Aj — ((/)/, A^'^-') for the fields, where the indices 
I, J, . . . label both scalar and vector fields; for scalar fields we use capital latin indices I,J,... and for vector fields 
we use lower case latin indices z, ... for vector components and another set of supra indices a, &, . . . to label the 
vector field |3] . In the following example we consider models involving only one vector field so, we can drop the supra 
index a, &, . . . For the derivatives of the amount of expansion N, we use: 



ON 



Nij - 



dAjdAj 



etc. 



(23) 



For simplicity, we shall consider models in which there is only one vector field so there is only one preferred direction 
given by the unit vector = A^i/(7VfeiVfc)^/^; however, the generalization for several scalar and vector fields should 
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be straightforward (for more details, see Ref. flT'). We will also impose the following conditions over the derivatives 
of TVH 



Nj oc Ay, Njj oc 5 



IJ- 



(24) 



The second condition impose that the interactions between scalar and vector fields are null or at least are negligible 
compared to the scalar-scalar and vector-vector interactions. Some of the most popular inflationary models in the 
presence of vector fields such as vector curvaton [52]-|Mji vector inflation [48J, and the hybrid inflation with coupled 
vector and scalar "waterfall field" where C is generated at the end of inflation [5T] are indeed parametrized as we 
present here. These cases certainly obey the conditions stated in Eq. (24). With the notation and conditions described 
above, the power spectrum of C is parametrized as follows [49j : 



Pc(k)-Pr(fc) (l+5c(k-fi) 



(25) 



where 



P^°{k)^NaNbPab{k) 



(26) 



is the isotropic part of the power spectrum with Pab being the power spectra of the field perturbations. The dimension- 
less power spectrum is related to the spectra of scalar and vector perturbations through 2TT^Pab{k) = k^^'^^"'''' Pabik) ■ 
For simplicity we assume that all the perturbations have approximately the same spectral index so that we take 
nab ~ ns- In any case, the results we exhibit here are not affected by the inclusion of this index and we could 
neglect it. The vector k is the unit vector along the wavevector k, the parameter g,^ quantifies the level of statistical 
anisotropy in C, and the vector fi specifies the direction of this anisotropy. In general, the statistical anisotropy 
parameter should acquire some level of scale dependence; however, in order to simplify our analysis, we consider 
here a scale-invariant anisotropy parameter. 

The bispectrum, expanded in g^, is parametrized as follows]^ 

i3c(kl, k2, kg) = (1 + ei)Sr(fcl, fc2, h) 



1 + ffc/ 



E (fc;fc„0-(3-"=) (n.k(,))2 + (n.k(„,))2 

l<7n 

E (fciA;„)-(3-"=) 



-2 p Km 

9C J 



E (fc/fc„)-(3-"=)(n • k(,))(n • V))(k(,) • k, 



(m)J 



l<.7n 



B'^iki, k2, fcg), 



(27) 



where k(;)j, with i = 1, 2, 3, denote the components of the vectors ki, k2, ka. 



5C = (ttif ) 5C, and / 



NiNiV 



+ 



(28) 



jy jy^y, — 1 . lu thc ab o vc expressions, V+ is the parity-conserving power spectrum 

defined as V+ = ^^^^^ where Vr^l are the right and left vector field spectra in a circular polarization basis (for 
details, see Refs. [35]). 



^ These conditions are rather typical when the background space undergoes isotropic expansion 1331 . One could argue, in principle, that 
these conditions still hold even when the background undergoes anisotropic expansion given that the potentials of the fields obey suitable 
slow- roll conditions. 

^ See Ref. 50 for a detailed discussion of the type of models that can be parametrized in this way. In that reference, the authors consider 
a scale-dependent anisotropy parameter 5^(fe) and even present a more general parametrization which involves non-Abelian gauge fields. 
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Finally, the trispectrum, expanded in qq, is parametrized as follows: 

r^"-^(ki, k2, k3, k4) = (1 +e2)Tf°(fci, fc2, fcg, k^) 



E [kikM-'^^-^'^^ (n.k(,))2 + (n.k(„))2 + (n.k(;,))^ 

l<.7n,s^m.l 

E (fcifc™fcz,)-(3-" = ) 



E {hkmkis) (n-k(i))(n-k(„,))(k(,) •k(„)) 

9c.h 



E (fci/c„fc;s)~^^""=)(n • k(,s)) (n • k(„))(k(„j) • k(,^)) + (n • k(i))(k(;) • k(;s)) 
9c.h 



9c h ' 



E {kikmkisy^^^^'^Hn ■ k(,))(n • k(„))(k(i) • k(,^))(k(,„) • k(,,))' 



E (fczfc„,fc,,)-(3-"=) 



(29) 



where we use the notation kmn — I km + k„|, h 



and 



fc2, A:3, fc4) = NaN^N^aN^fVacVbeVdf ^^kM''"''""'^ 

l<im,Sy^m,l 



(30) 



B. The non-gaussianity parameters 



The non-Gaussianity parameters /nl and tnl defined through Eqs. Q and ^ can be obtained directly from Eqs. 
( 25 1 , ( 27 ) , and ( 29 ) . Since our intention is to describe cases in which the fourth variety of the S Y consistency relation in 



Eq. (21 1 is violated, a simple minded example is enough for our purposes; then, we consider a single vector field so that 
gc = Qc and f — h = 1. For our evaluations, we use the parametrization described in Fig. [ij As shown in the figure, the 
wavevectors k^ are restricted to form a warped polygon due to the statistical homogeneity and the vector n is described 
by its polar and azimuthal angles 9 and respectively. In sum, the bispectrum of C and the level of non-gaussianity /nl 
are functions of the lengths of the triangle formed by the wavevectors, the angles that parametrize the orientation of n, 
and the level of statistical anisotropy g^; this is: /nl = /nlI^i, ^2,^3, 0, (j), g^)- Analogously, the trispectrum of and 
the level of non-gaussianity tnl have the following functional dependence: tnl = ''nl(^1j ^2: ^Sj ^127 ^14: ^, 0, 5c)- 
Thus, following Fig. [l] /nl = f-NL{k,x,y,9,(t),gc) and tnl = Tyi-L{k,x,y,p,q, z,0, cj), g^). 



Results 



In order to test the fourth and fifth varieties of the SY consistency relation, Eqs. (21 1 and (22 1, we have to take the 
squeezed limit of /nl, Eq- ([8|, and the collapsed limit of tnl, Eq. ([6]). The squeezed limit of /nl has to be evaluated 
both when ki — >■ and when ks — > 0. Then, looking at Fig. [ija) and using the parametrization f^hik, x, y, 9, (/), g^) 
we take the first limit as: 



/NL(ki) = lim /NL(ki,k2,k3) = /NL(ki, -ki, 0) = fj^Lik,k,0,9,(f>,gc) 



(31) 



For the second limit, we look at Fig. [ijb) and use the parametrization f]<si^{k,x,y,p,q, z,9,(l>, g^): 
/NL(k3) = lim /NL(ki2,k3,k4) = /NL(0,k3, -k3) = f^i^{k,k,0,p,q,z,9,(l),gc) . 

ki2->-0 



(32) 



As a result, both limits are independent of the scale k. On the other hand, in the collapsed limit of tnl, when ki2 — >■ 0, 
we obtain: 



rNl(ki,k3)= lim TNL(ki, k2, k3, k4) = TNL(ki, -ki, ks, -k3) = TNL(fc, fc, 0,p, g, z, 6*, 5<;) , (33) 

ki2->0 



9 



=(sert 8coscp , sen^sen cp , cos 9} 



k,=k 

k^ = i(x-kf+y~ 



k, = {x-k,-y,0) 



k,=(k,0,0) 



k^ = {-x,y,0) 



= { sen 6 cos cp , .?e/7 9 sen cp , cos 9 




(a) Parametrization of the bispectrum of and /nl- (b) Parametrization of the trispectrum of and tnl- 

FIG. 1: The figures (a) and (b) represent the parameters employed to calculate the correlation functions of ( and 
the non-gaussianity parameters. Aside of the wavevectors, when statistical anisotropy is taken into account, one 
needs a unitary direction of a preferred direction. Here, this direction is represented by the vector n which is 

parametrized using spherical coordinates (^, (/>). 



which is also /c-independent. Finally, the collapsed limit evaluated at ki is obtained from the previous expression in 
the following way: 



TNL'(ki,ki)= lim TNL(ki,-ki,ki,-ki) = rNL(fc,fc, 0,^,0,0,6*, 0, 



9C) 



(34) 



which is again fc-independent. The results of (kO/^'f (k3))/r™"(ki, kg) and (|/^l^(ki))VT-"(ki, ki) evalu- 

ated in a particular configuration are shown in Fig. [2j The results shown in the figure correspond to the case in which 

36 /'i'(kl)^ 



36/^l'(kl)/;«J(k3) 



25 r 



collapsed . 




25 T'°"(kl, kl) 

2.5 : 

2.0 
1.5 
1.0 
0.5 



-0.5 ; 
-1.0 ; 



-0.5 



kl=k 

<p=7r/3 



8( 



(a) There exist configurations in which the fourth variety of (b) The fifth variety of the SY consistency relation (Eq. [22| is 

the SY consistency relation (Eq. [2l]) is violated. Here, as an preserved in any wavevector configuration. As an example, 
example, the configuration p = q = z^x^k,e = it/6, and this figure shows the evaluation of Eq. ([22| in the 

= 7r/3 is showed. configuration ki = k, 9 = n/6, and <j) = n/3. 

FIG. 2: The fourth and fifth varieties of the SY consistency relation under test. 



all the scales are equal, this is: p = q = z = x = k, and the unit vector along the preferred direction is characterized 
by = 7r/6 and (j) = 7r/3. Both quotients in this configuration are independent of the scale k but exhibit a non trivial 
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dependence on the orientation of the unit vector n and the level of statistical anisotropy g^. In Fig. [2]ja) , we have an 
explicit example of the violation of the fourth variety of the SY consistency relation, Eq. (21 1, while in Fig. |2jb), we 
have an example of the preservation of the fifth variety of the SY consistency relation, Eq. (|22p. The case represented 
in Fig. [2]ja) is the relevant result for our purposes because this is the only proof that we need in order to show that the 
inequality (21 ) is not a generic feature of models with statistical anisotropy. There are many other (actually infinite) 
configurations in which violations of the fourth variety, Eq. (21 ), are obtained. The result shown in Fig. [2jb), on the 
other hand, is only a particular case of the fifth variety, Eq. (22 1, but we can see that this relation is respected in any 
configuration of the wavevectors and the vector fi and for any value of the level of statistical anisotropy . 

It is worth noticing that we do not consider just the values of in agreement with current observational bounds 
{g(; = 0.290 ±0.031) [lH| since it is not our intention to restrict the analysis to realistic models of statistical anisotropy, 
being that the reason why in our evaluations we indeed consider larger values of 5^, from to 10; this is done because 
it helps us to realize that no matter how large the statistical anisotropy in ( can be, the fifth variety of the SY 
consistency relation, Eq. ( |22[ ), is always respected, which means that this relation is independent on the assumption 
of statistical isotropy in the primordial curvature perturbation. On the other hand, it is also worth noticing that, 
even for values of in agreement with current observations, we can obtain violations of the fourth variety of the SY 
consistency relation, Eq. (21), within percent level size. 



V. THE PRIMORDIAL CURVATURE PERTURBATION GENERATED BY SCALAR AND VECTOR 
FIELDS: THE SIXTH VARIETY OF THE SY CONSISTENCY RELATION 



We briefly present in this section the latest variety of the SY consistency relation which was formally demonstrated 
in a quite recent paper by the authors of this work |47j : 



'NL 



/•ISO 

./NL 



(35) 



In the previous expression, t^£ and /j|£ correspond to the isotropic pieces (i.e., they do not depend at all on the 
wavevectors) of the levels of non-gaussianity in models where the primordial curvature perturbation is generated by 
scalar and vector fields. One could think that t^£ and f^^, being isotropic pieces, have contributions only from the 
scalar fields, but, as demonstrated in [37], they receive contributions from the vector fields as well. 

As done in Section 1, we will enumerate the conditions under which the sixth variety of the SY consistency relation 
is satisfied: 

• Condition 1: The calculation of /jlf£ and r^£ is performed at tree level in the diagrammatic approach of the 5N 
formalism [H [5] . 

• Condition 2: The expansion is assumed to be isotropic so that the field perturbations (actually the scalar 
perturbations that multiply the respective polarization vectors (if any) in a polarization mode expansion) are 
statistically isotropic. 

• Condition 3: The infiationary dynamics is driven by any number of slowly-rolling vector and scalar fields. 

• Condition 4'- The fields involved are gaussian. 

• Condition 5: The field perturbations (with the proviso expressed in Condition 2) are scale-invariant. 

As discussed previously, /nl and tnl are not dir ectl y comparable in the general case because they are functions of 
the wavevectors; however, the expression in Eq. (35) is valid since, under the conditions just stated above, /jlf£ and 
T^£ are scale- invariant [55] . 



VI. CONCLUSIONS: THE SY CONSISTENCY RELATION(S), NOT JUST ONE BUT SIX 

When first introduced back in 2008 [S] , the SY consistency relation was considered as a useful tool to discriminate 
among different classes of models for the generation of the primordial curvature perturbation; however, the quite strong 
conditions under which this relation was proven led to several authors to wonder whether such conditions could be 
relaxed in the search of a more accurate and clean way to employ this consistency relation as a discriminator tool. 
Since that time, six different varieties of the SY consistency relation have been implicitly introduced in the literature. 
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one of them indeed introduced in the present paper. In this work, we have collected the six varieties and showed them 
explicitly as well as the conditions under which they are valid. The first and the sixth varieties, Eqs. ([!]) and (35 1, 
are quite similar in the sense that the non-gaussianity parameters /nl ^^nd tnl are scale-invariant; they are different 
in the sense that the first variety involves the whole /nl and tnl while the sixth variety involves only their isotropic 
pieces. The second and third varieties of the SY consistency relation, Eqs. (18) and (19), are very similar since they 
involve an integration over the momentum space but are different because the integration in the second variety is just 
in a narrow band of wavevectors whereas in the third variety the integration is over the entire momentum space while 
introducing the power spectrum of the primordial curvature perturbation as part of the integrand; anyway, because of 
the integrations, it looks quite difficult to employ these varieties of the SY consistency relation in order to discriminate 
among different classes of models. What is interesting however of these varieties is that they are valid even if there 
is neither statistical isotropy nor scale-invariance: that is why the fourth variety in Eq. (21) is so interesting, being 
just the best-defined generalization of the first variety and whose form is easily inspired from the second and third 
varieties; it is for sure valid for models involving only slowly-rolling Lorentz-invariant scalar fields, as long as the proof 
by Kehagias and Riotto [M] turns out to be valid at all orders in perturbation theory, but for sure fails in some model 
that involves non-trivial degrees of freedom; indeed, we have shown in Section [IV] the first counterexample: the fourth 
variety is strongly violated in a quite generic class of models that involves vector fields and that include the vector 
curvaton model proposed by Dimopoulos [22]- [24], the vector inflation model proposed by Golovenev, Mukhanov, and 
Vanchurin [48 , and the hybrid inflation model with coupled vector and scalar "waterfall fleld" where C is generated at 
the end of inflation proposed by Yokoyama and Soda [21 . The latter results sounds ubiquitous for any variety of the 
SY consistency relation, that does not involve integrations over momentum space, when there is no scale-invariance; 
nonetheless, we have been able to introduce the fifth variety of the SY consistency relation in Eq. (22) whose validity 
is only restricted to the existence of statistical homogeneity and when /NL(ki, k2, ks) and rNL(ki, k2, ka, k4) are 
calculated in the squeezed and collapsed limits respectively; indeed, this variety stands the test of the class of models 
described above as shown in Section |IV| perhaps the strongest conclusion we can extract from all this analysis is a 
simple question: what will we do if the fifth variety of the SY consistency relation is observationally violated? It 
would imply statistical inhomogeneity, which in turn implies that, although we dispose of very good data and very 
good theory, we cannot compare them to each other (4i il2] : this situation would be a real shake in the foundations 
of cosmology as a science [i^ . 
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